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Abstract 

Programs that calculate observables in quantum chromodynamics at next-to-leading order typi- 
cally generate events that consist of partons rather than hadrons - and just a few partons at that. 
These programs would be much more useful if the few partons were turned into parton showers, 
which could be given to one of the Monte Carlo event generators to produce hadron showers. In 
a previous paper, we have seen how to generate parton showers related to the final state collinear 
singularities of the perturbative calculation for the example of e + + e~ — >• 3 jets. This paper 
discusses the treatment of the soft singularities. 
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I. INTRODUCTION 



This is the second of two papers concerning how to modify a next-to-leading order (NLO) 
calculation in quantum chromodynamics by adding parton showers in such a way that when 
the generated events are used to calculate an infrared safe observable, the observable is 
calculated correctly at next-to-leading order. In the companion paper [1], Kramer and the 
present author explain the aspects of this problem that relate to the divergences that appear 
in perturbation theory when two massless partons become collinear. This part of the method 
is quite simple and easy to understand. The treatment of the divergences that appear when 
a gluon becomes soft is a bit more technical and will be given in this paper. 

We find that it it is conceptually simplest to match showers to a NLO calculation if 
the NLO calculation is performed in the Coulomb gauge. That is because the collinear 
divergences of the theory in a physical gauge like the Coulomb gauge are isolated in the 
cut and virtual self-energy diagrams - just the diagrams that you would draw to illustrate 
a parton splitting to two partons. Accordingly, the present research program was initiated 
in Ref. [2], which shows how to do NLO calculations in the Coulomb gauge. I should point 
out that, although there is a certain conceptual simplicity that emerges through the use of 
a physical gauge, one can use any gauge, as the recent work of Frixione, Nason, and Webber 
on the same subject demonstrates [3, 4]. 

I begin with a brief restatement of the problem. Consider the calculation of an exper- 
imental observable that describes a hard scattering process involving strongly interacting 
particles and suppose that the chosen observable is infrared safe, which makes a perturba- 
tive calculation for the observable sensible. The simplest sort of calculation for this purpose 
consists of calculating the hard process cross section at the lowest order in a s , call it af , at 
which it occurs. However, one often performs a next-to-leading order calculation, including 
terms proportional to af and af +1 . Then the estimated error from yet higher order terms 
(which are usually unknown) is typically smaller than with a leading order calculation. 

The NLO calculations are typically in the form of computer programs that act as Monte 
Carlo generators of partonic events. They produce simulated events with final states / 
consisting of a few partons with specified momenta. For example, for three-jet production 
in electron-positron annihilation, there are three or four partons in the final state. Each 
event comes with a weight w n . Then the predicted value for an observable described by a 
measurement function S(f) is 

1 N 

n=l 

A serious weakness of most current NLO Monte Carlo generators is that they produce 
simulated final states that are not close to physical final states. This contrasts sharply with 
the highly successful leading order Monte Carlo event generators such as Pythia, Herwig, 
and Ariadne [5-7], which do generate realistic final states. What could be done to generate 
realistic final states in an NLO calculation? One might think of using the three and four 
parton final states in the NLO partonic generators to generate a parton shower from each 
outgoing parton. This gives a final state with lots of partons. Then one could use one of the 
leading order Monte Carlo event generators to turn the many partons into hadrons. There 
should not be a serious problem with the hadronization stage, since this is modelled as a long 
distance process that leaves infrared safe observables largely unchanged. The problem lies 
with the parton showers. Here a high energy parton splits into two daughter partons, which 
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each split into two more partons. As this process continues, the virtualities of successive 
pairs of daughter partons gets smaller and smaller, representing splittings that happen at 
larger and larger distance scales. The late splittings leave infrared safe observables largely 
unchanged. However, the first splittings in a parton shower can involve large virtualities. 
They represent a mixture of long distance and short distance physics. Thus one must be 
careful that the showering does not reproduce some piece of short distance physics that was 
already included in the NLO calculation. To be a little more precise, if one expands the 
prediction of the calculation including showers in a power series in a s (Q), where Q is the 
hard process scale, in the form 

l = C af + C ia f +1 + C 3 af +2 + ---, (2) 

then our goal is to insure that the coefficients Cq and C\ are exactly the same in the 
calculation with showers as in the purely perturbative NLO calculation. 

The algorithms that we propose apply to showers from final state massless partons. They 
are illustrated using the process e + + e~ — > 3 jets in quantum chromodynamics. 1 While our 
research program was underway, Frixione and Webber succeeded in matching showers to a 
NLO calculation for a problem with massless partons in the initial state (but no observed 
colored partons in the final state) [3]. Subsequently, Frixione, Webber, and Nason have 
extended this method to a problem with massive colored partons in the final state [4]. Their 
method is similar to the method presented in Ref. [1], or more precisely to a variant of 
that method in which the first level of parton splitting, based on Coulomb gauge splitting 
functions in Ref. [1], is instead based on the splitting functions of Herwig. The soft gluon 
effects treated in this paper are not part of the algorithms of [3, 4]. Instead, these references 
employ what amounts to a soft gluon cutoff. 2 There has been other work on this problem 
[8]. That of Collins [9] is particularly instructive. 

II. THE PROBLEM OF THE SOFT SINGULARITIES 

In this section, I describe the particular problem left over from Ref. [1] that remains to 
be addressed in this paper, recalling just enough of the argument of Ref. [1] to set up the 
issue. 

Consider one of the Born graphs for e + + e~ — > 3 jets. Three partons emerge into the 
final state. The idea presented in Ref. [1] is to let each of these partons split into a pair of 
partons, as illustrated in Fig. 1. Two features of the splitting are important. First, when the 
virtuality q 2 of the pair is small, the splitting probability is proportional to (P(x)/q 2 )dq 2 dx 
where x is the share of the momentum carried by one of the daughter partons and P(x) is the 
appropriate Altarelli-Parisi splitting function. Second, the collinear singularity at q 2 — > 
is damped by a Sudakov factor with the behavior exp(— a s c log 2 (g 2 )) for q 2 -> 0. Thus 
each of the three partons makes a jet, but in the limit of small a s , each jet is usually very 
narrow and appears in an infrared safe measurement like a single massless parton. There 
is a qualifying adverb "usually" here. A fraction a s of the time, one of the splittings has a 



1 I discuss e + + e — » 3 jets but not e + + e — » 2 jets. Thus the present calculation is to be applied with a 
measurement function that give zero for events that are close to the two-jet configuration. 

2 The parameter C\ above is then modified by an amount proportional to a positive power of the cutoff 
parameter (3, so as long as (3 is chosen small enough, the cutoff effects in [3, 4] are negligible. 
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FIG. 1: Primary parton splitting from Ref. [1]. The filled circles represent graphs for the Born am- 
plitude and complex conjugate amplitude. Each of the partons emerging from the Born amplitude 
splits into two partons with a vertex, represented by the squares, that includes a Sudakov suppres- 
sion factor. Each of the six daughter partons undergoes further, secondary, splittings and enters 
the final state as a complete shower. The secondary splittings are represented by the diamonds. 

substantial virtuality and we get a four-jet final state. Thus there is an order af +1 effect in 
which some probability is removed from the three-jet final state and given to a four-jet final 
state. In a purely NLO calculation, this effect is included in the af +l graphs. Thus, to keep 
the calculation correct to NLO, we subtract these probabilities from the af +1 perturbative 
graphs. 

The subtractions are required for NLO consistency of the calculation, but they have 
another important effect. Since the splitting functions have just the structure found in the 
order af +1 graphs in the collinear limit, the subtractions cancel the collinear singularities 
(and the collinear x soft singularities) of the order af +1 graphs. The real and virtual 
singularities are separately cancelled, whereas without the subtractions the real and virtual 
singularities are left to cancel against each other (which, of course, they do 3 ). If there were 
no other singularities, we could let the three and four parton final states of the order af +1 
graphs develop into parton showers, which would affect the result of the calculation only at 
order af +2 . The calculation would be numerically stable because the integrals for the three 
and four parton states would be separately finite. 

But, of course, there are other singularities, namely the singularities from emission of soft 
gluons at wide angles from the three hard partons, as well as the corresponding singularities 
in the af +1 virtual graphs. These singularities cancel between the order af +1 graphs with 



3 We employ a calculation [10-12] in which the integrals over the three-momenta in virtual loops are 
performed numerically. Then this cancellation happens automatically point by point in the integration. 
In other methods, the virtual integrals are performed entirely analytically. In this case, the cancellation 
mechanism is less simple. 
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three and four final state partons, but we cannot let these two kinds of final states develop 
separately into showers without losing this cancellation and thus encountering numerical 
instability. 

This is not just a technical problem to be somehow sidestepped. The emission of a soft 
gluon at order af +1 is a real physical effect. It is an effect that does not fit a picture of 
separate evolution of the three jets because the soft gluon is sensitive to the whole final 
state. That is, the probability to emit a soft gluon does not break up into the sum of three 
independent probabilities, one for emission from each jet. 

There is a natural way to deal with this. We can consider the three jets to form an antenna 
that radiates a soft gluon. The radiation probability is proportional to (f ($,(/>)/ E)dEdQ, 
where E is the energy of the gluon, 9,4> are its angles, and f(9,<f)) is a function that is 
determined by the directions of the three jets. We can account for the normalization of the 
radiation probability provided by the virtual graphs by providing a suppression factor with 
the behavior exp(— a s c \ \og(E)\) for E — > 0. This is depicted in Fig. 2. Thus a soft gluon is 
always emitted, but when a s is small the gluon is usually too soft to matter to an infrared 
safe measurement. However, a fraction a s of the time the gluon can have a substantial energy 
and we get a four-jet final state. Thus there is an order af +1 effect in which some probability 
is removed from the three-jet final state and given to a four-jet final state. In the purely 
NLO calculation, this effect is included in the af +1 graphs. Thus, to keep the calculation 
correct to NLO, we subtract these probabilities from the af +l perturbative graphs. 4 Of 
course, the required subtraction terms will be just the soft singularity subtractions that we 
needed to stabilize the af +1 contributions if we add showers separately to af +1 graphs with 
three and four-jet final states. 

Does the necessity to add soft radiation from a three-jet antenna mean that the standard 
parton shower Monte Carlo programs such as Herwig, Pythia, and Ariadne need to be 
rewritten in order to be compatible with an NLO calculation? Not at all. We do not need 
soft radiation with its correct angular dependence at every stage of a parton shower. All that 
we need is to radiate one soft gluon from the hard event with its correct angular dependence. 
This radiation, along with the initial splittings of the three hard partons, can be considered 
to be part of the NLO calculation, so that the Monte Carlo programs do not need any 
modification at all. I will add more comments on this subject in the conclusions section of 
this paper. 

III. SOFT GLUON EMISSION 

In Ref. [1], we saw that one can turn "Born x [1 + real — \virtual\}" into "Born x 
real exp(— \virtual\)" with respect to self-energy graphs, which contain the collinear and 
collinear x soft singularities in the Coulomb gauge. In this section, we begin a detailed 
description of the singularities generated by the emission of soft gluons at wide angles to the 
outgoing partons. Up to now in our description, these singularities cancel between real and 
virtual emissions in the form "Born x [1 + real — \ virtual\}." In the following sections, we turn 
this into the form "Born x real exp(— \virtual\y together with some left over nonsingular 
terms. 



As we shall see, some of the adjustment is to be made in the wide angle part of the splitting functions 
used to generate the splittings of the three hard partons. 
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FIG. 2: Primary parton splitting with soft gluon radiation added. The graphical symbols of Fig. 1 
are also used here. The extra gluon coming from the circles on the right and left represents the soft 
gluon radiated from the three jets. As with the other partons, the soft gluon undergoes secondary 
splittings and generates a shower, represented by the diamonds. 

Consider a cut Born graph in which the final state has a quark (momentum qi), a gluon 
(momentum q%) and an antiquark (momentum q%). These partons will make collinear show- 
ers, but we ignore that for now. What concerns us here is that these outgoing partons act 
as an antenna that radiates soft gluons. Let a soft gluon with momentum I be emitted by 
parton i and, in the complex conjugate amplitude, absorbed by parton j. 

We need to specify the kinematics in a definite way. We choose to route the momenta 
for the {i,j} term so that the momenta {px^\p2 > P3 } carried by, respectively, the quark, 
the gluon, and the antiquark in the final state are given by 

= qk-\{^ + h 3 % (3) 

Thus for a graph with i = j, the momentum k before the soft gluon emission is k = q and 
afterwards it is p = q — I. For a graph with % 7^ j, the momentum k before the soft gluon 
emission is k = q + ~ 1 and afterwards it is p = q — 1 1. 

A. The soft gluon approximation 

We now specify the approximation to be applied for the soft gluon emission [13]. Let a 
parton with momentum k emit a soft gluon with momentum I and enter the final state with 
momentum p. Then P = —p^. The momenta p and k are approximately in the direction 
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— * 

of a light-like vector u. We take u = (1, q) where q = q/\q\, and q is \(k + p) in the case of 

— * 

an % ^ j graph or k in the case of an i = j graph, as specified in Eq. (3). 

The approximation is based on the fact that the soft gluon sees a current J M that is 
approximately proportional to w M . Thus if we denote the polarization vector of the gluon by 
e M , we can replace 

J-e^J-l^-. (4) 
u ■ I 

This replacement becomes exact in the limit in which J 11 approaches Aw M for some constant 
A. It is important the the momentum / not point along u. That is, the gluon momentum 
should be small so that both k and p and thus J are approximately proportional to u, but 
the gluon momentum should not be collinear with the momentum of the emitting parton. 

For the graph in which the soft gluon goes from parton i to parton j, the factors u ■ e/u ■ I, 
together with appropriate color matrices, give a factor (in the Coulomb gauge) 



u 



13 {u?Q{u%) ~ P (1 - 1 ■ ^ (1 - 1 ■ qj) ' 
Here D^ V {1) is the Coulomb gauge numerator factor 

d» v {i) = + 1 [-in v - ln v + m v \ , (6) 

where / is (0, 1 1 , 1 2 , 1 3 ). The constants are the color factors for this process, 

N c 

O12 — O21 — O23 — O32 — — , 

1 



C\z — C31 



2iV 
C11 = C33 — Cf, 

C22 = C A . (7) 
What about the factor J -P. In the case of emission from a quark, we have (since p 2 = 0) 

(8) 

so we simply include a factor j> for the final state quark. The final state phase space associated 
with the original graph has a factor l/(2|p|). We retain this factor without change. 

For emission from an antiquark we have [noting that the momenta p and q are directed 
against the fermion number flow and that we should extract a factor (—1) to associate with 
the color factor for the graph] 

(-l)x fc£(f-fl(-fl = -ji, (9) 

so we simply include a factor — j> for the final state antiquark. The final state phase space 
associated with the original graph has a factor l/(2|p|). We retain this factor without 
change. 
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The case of emission from a gluon is a bit more complicated. We denote the tensor 
associated with the vertex function by 



V a ^(k A , k B , he) = g a \k\ - kl) + cyclic permutations. (10) 

Then the three gluon vertex function with momenta {Au, &b, kc} directed into the vertex is 
(—ig)(—if a b c )V al3l (kA, k B , kc)- The vertex dotted into (k—p) times the adjoining propagator 
and the propagator numerator for the final state gluon, with the —ig and the color matrix 
-ifabc removed, is 

D(kY 

d{pX [vf(-p,k, P -k)(w-pi)] 

DikY 

= -D{ P y a [{k 2 g a P - k a k p ) - ( P 2 g a ? - P a P y 

Dm 



,2 



= -D( P y a (k 2 g a P -k a k p ) 
k z 

= D(py a (g- + ±-(p« + hkn 
k 2 

= D(py a (9 au + Y 2 i a kn 

-> D( P r. (ii) 

Here in the second to last step we have noted that D(p)£p a = 0. Then in the last step we 
drop the term proportional to I so as to simplify our soft gluon approximation. Thus the net 
result of the approximation is to include, in addition to the factor u • e/u • Z, a factor D{pY u 
for the final state gluon. The final state phase space associated with the original graph has 
a factor l/(2|p|). We retain this factor without change. 

It is important that the factors in Eqs. (8), (9) and (11) revert to the factors in the Born 
diagram when I — > 0. 



B. The added terms 

We can now apply these ideas. We write the contribution to the cross section from a 
given cut Born level graph as 

J[Born] = jdqi J dq 2 J dq 3 5{^q^j G 3 (q 1 ,q 2 ,q 3 ) 

x/s(5i, 52, 9s) £3(91,92,93) ■ ( 12 ) 

Here S 3 is the final state measurement function [see Eq. (4) of Ref. [1]] and f 3 contains 
the factors from the Feynman rules corresponding to the three final state particles, which, 
suppressing the Dirac and Lorentz indices, are 

fin ft n \ — ^ D(g2) ~^ m\ 
A \Qi\ *\Q2\ *\Q3\ 
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The function G3 contains everything else, including an amputated Green function times a 
hermitian conjugate amputated Green function. With this notation, we write the terms to 
be added for real soft gluon emission as 



T[soft, real] = Jdq\ J dq 2 J dq z 5{^q^ G z {qi,q 2 ,q\) 

9s f dl _ 
(2tt) 3 J 2\l \ I 



n 2 r dl 1 

x jtf iw\ P 9{P K Ms ° ft) E F ^ ^ ^ 

x fM ij \ pP) pE pE I) , (14) 



where are defined in Eq. (3) and S4 is the measurement function for four final state 
particles. 

The factor did < M s 2 oft ) is inserted to suppress contributions from momenta / that are 
outside the range of validity of the soft gluon approximation. We take 

M soft = A softv /i^(l - t ), (15) 

where y^so and t are, respectively, the cm. energy and the thrust of the final state with 
parton momenta {q\, q 2 , cfe} and where we take the parameter A so f t to be 1/3. 
The functions are, for i ^ j, 

Fij = C l3 * A" jjMjM, (16) 
as specified in Eq. (5). For % = j we begin with a factor 

c ,l4i!! = c ,li!i, (17) 

"(1 -;•&)' l-l-q. 
We use Eq. (7) to break Cu up into 

3 

Cu = — E Cia- (18) 

a=l 

Then in the term proportional to Cj„, we impose the condition that the angle between I and 
qi should be greater than the angle between q a and q}. That is, our soft gluon for % — j 
should be soft but not collinear, and we impose a cut to enforce that condition. Thus we 
take 
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l-l-qi 



{-!) C ia 6(l ■q i <qa-qi 



(19) 



It will be useful to have at hand an abbreviated notation based on the definitions 



Ro(qi, qs) = G 3 (qi, $2, qs) fs{qi, <h, qs) 5 3 (gi, <f 2 , qs) (20) 

and 

Rv$ = Gsi^^qs) HpEpEpE) SM ij \p? j \pEl)- (21) 
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With this notation, the Born and soft-real contributions are 

J[Born] = Jdqi J dq 2 J dq 3 S(j2 Qi) Ro(Qi, 92, «0 (22) 



and 



^ S o ft d\i\ r d 2 i 



X[soft, real] = J 'dq\ J dq 2 J dq 3 5 (J2 Qi) J Q -=) J 



l\ J 4tt 

X E~ F ij( 1 '' 9i, 92, qa) Rijih 9i, 92, 9s)- (23) 
ij n 

The indicated integration over I is an integration over the unit sphere. 

We define a corresponding virtual soft gluon contribution that will cancel the real soft 
gluon contribution for small / as long as the measurement function is infrared safe, so that 

— * 

Rij — > R for I — > 0: 



X[soft, virtual] = — J dqi J dq 2 J dqs S^Yoi) J 



M *oit d\l\ f dH 



\l\ J 4?r 



X Yl — 9i> 92, 9s) Ro(Qi, 92, Qs)- (24) 

— * 

Conveniently, the integral over the angles of / in X [soft, virtual] is simple: 

f d 2 l - 2 3 

/ 7- 2 ?i,?2,93) = E E (-Cij)^ ~ Qi ■ Qj) 
47r ij i=lj=i+l 

N r 1 
= -^[(1 " 9i " 92) + (1 - 93 • 92)] - 2^(1 - 9i " 93). (25) 

The integral is largest when q\ ■ q% = 1, which implies q\ ■ q 2 — <?3 • q 2 — — 1- Then the integral 
is 2iVc = 2Ca- The integral is smallest when q± - q 2 = 1, or else q^ ■ q 2 = 1. Then the integral 
is N c - l/N c = 2C F . 

We will, in due course, see how to exponentiate the added contributions X[soft,real] + 
X[soft, virtual]. First, however, we must subtract these terms from X so that the cross section 
remains unchanged. 



IV. ADJUSTING THE COLLINEAR TERMS 

The terms in Eqs. (23) with i — j represent soft gluon approximations to the cut self- 
energy diagrams for outgoing partons i. Thus we should subtract them from the exact cut 
self-energy diagrams. Consider, for example, % — j — 1, corresponding to a self-energy 
insertion on the outgoing quark propagator. This contribution has the form given in Eq. (7) 
of Ref. [1], 

^'l = IW\ Mf? iV'f ^M 3l ,(f, x A)R(f, X A)}, (26) 
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where q = q± and the notation, including the variables (q 2 ,x,<p), is explained in Sec. Ill of 
Ref. [1]. We can translate the {1,1} term in Eq. (23) into the notation used in this equation 
by using the change of variable 



/ 



dl 



2\q-l\2\l 



| rOC rl pTT ] 

= mJ° df L dx L d(l) TTA 



where 
Then 



1 + A = y/l + f/q 2 . 



(27) 
(28) 



Ill[soft>real] . | * ^ e(P K <ft 



'dq 



«2 ,1 



*" d0 1 g 2 a s 



x- 



1 + 

1 - Z • g 



:-l)X)Ci o 0(f. g<g„-g 



Now we can use 



a = l 



f 2 <? 2 



-l)R{q\xA) ■ (29) 



Z 2 = — (A + 2x) 2 
4 

g 2 = g 2 A(2 + A) 
f , 2x-A + 2xA 
l ' q= A + 2x 



(30) 



to rewrite this as 
Xn[soft, real] 



dq 

W\ 



Tr< 



•2 -i 



> <ig 



o g^ jo 



7T27T 1 + A 



fl(A + 2a;<2M soft /|g|) 



x 



2x - A + 2rc A 



x- 



yx 



A + 2x 
1 + A/2) 2 



<Qa-q 



- /) R(<i 2 ,x,> 



2tt (A + 2:r) 2 1-x 
Thus we can subtract J n [soft, real] by replacing M g / q (q 2 , x, 4>) by 



(31) 



(32) 



where 



:-Di;ay 2 *- A+2 * A 



x 



8x 1 (1 + A/2) 

1 -x (A + 2a;) 2 1 + A 



a = l 

2 



A + 2x 



<Qa-q 



(33) 



The effect of this is to remove the part of M. that is singular when \l\ — > at a fixed angle 
greater than the angle between the quark and the other two outgoing partons. Note that 

— * 

\l | — > at a fixed angle corresponds to A — > and x — > with A/x fixed. 
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The term in Eq. (23) for emission from the antiquark line, X 33 [soft, real] is given by 
Eq. (31) with (^ — /) replaced by + /) and the index 1 replaced by 3. Thus we can 
subtract X 33 [soft, real] by replacing A4 g /g(q 2 ,x,(p) in the analogue of Eq. (26) by 

M g/q (f, x, 0) - {-i + /) P 3 so V, x), (34) 

with the function P| oft ( q 2 , x) given in Eq. (33) with i — 3. 

Finally, the term in Eq. (23) for emission from the gluon line, X 22 [soft, real] is given by 

J 2\q I Jo q z Jo J-nzn 

x ^ ^°V, x)D^{q - l)R,Af, x, 0), (35) 

with the function P 2 oft (<? 2 , x) given in Eq. (33) with % — 2. Thus we can subtract X 22 [soft, real] 
by replacing J\4 g 1 g (q 2 , x, 0) in the analogue of Eq. (26) by 

M^ g {q\ x, 0) - D^(q - I) Vf\q\ x). (36) 

(There is also a g — > gg splitting function M?£j g , but there is no soft singularity for g qq 
so this function does not get modified.) 

The terms in Eq. (24) with i = j represent soft gluon approximations to the virtual 
self-energy diagrams for outgoing partons i. Thus we should subtract them from the exact 
virtual self-energy diagrams. Consider the virtual self-energy insertion on the outgoing quark 
propagator. This contribution has the form given in Eq. (12) of Ref. [1], 

XMual] = /Hi Tr {- f f j\ lx £g £ V 3/ ,if, X ) iBt }. (37) 

The functions V a /b are given in the Appendix of Ref. [1]. When translated into the notation 
used in this equation, the {1,1} term in Eq. (24) is 

I„ [soft, virtual] = J Tr {_ J~ f fafjt |p PV , x) ^} , (38) 

where p^oft j g d e fi nec i j n Eq_ (33). Thus we can subtract Xn [soft, virtual] by replacing 

V 9/q (q 2 ,x)bj 

V g/q {q 2 ,x)-Vf\q\x). (39) 

Similarly, we can subtract X 33 [soft, virtual] by replacing V g /g(q 2 ,x) = V g / g (q 2 ,x) in the ana- 
logue of Eq. (37) by 

V g/q {q\x)-Vf\q\x) (40) 
and we can subtract X 22 [soft, virtual] by replacing V g / g (q 2 ,x) in the analogue of Eq. (37) by 

V glg {q\x)-Vf\q\x). (41) 

After the collinear splitting terms have been adjusted to remove their wide angle soft 
parts, we turn " Born x [1 + real — \virtual\}" for the collinear splittings into 11 Born x 
real exp(— \virtual\) v as explained in Ref. [1]. The only difference is that the meaning of 
"rea/" and "virtual" here have been slightly changed. 
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V. ADJUSTING THE NLO GRAPHS 

It remains to subtract the terms X^ [soft, real] and Xij [soft, virtual] for i ^ j from X. 
Recall that we include in the calculation the cut order af +1 graphs that do not have a cut 
self-energy subdiagram or a virtual self-energy subdiagram with the immediately adjacent 
propagator cut. (We also include diagrams with cut self-energy subdiagrams, but with a 
theta function that requires the virtuality in the self-energy subdiagram to be large.) The 
included graphs have soft gluon divergences, which cancel between real and virtual soft 
gluon emissions. We will subtract the approximate soft gluon contributions from these 
exact contributions, so that the soft gluon divergences are cancelled separately for the real 
and virtual contributions. 

A. Real soft gluons 

This is simple for the real emission terms, X[soft, real]. Consider a cut order af +1 graph 
with four parton lines cut. The cut lines are a quark and an antiquark plus either two gluons 
or another quark and antiquark. In the case of a quark and an antiquark plus another quark 
and antiquark, we do nothing. In the case of a quark and an antiquark plus two gluons, 
we designate one of the gluons as potentially soft. Then this cut graph corresponds to the 
contribution from one of the possible Born graphs to one of the terms in X[soft, real], which 
has two gluons in the final state, one of them "hard" and one soft. We simply subtract this 
contribution to X[soft,real] from the cut graph in question. Now we designate the other 
final state gluon as potentially soft, find the corresponding contribution to X[soft,real] and 
subtract it. 

This procedure has two effects. First, we use up all of the terms in X[soft,real] that we 
needed to subtract from X. Second, the cut next-to-leading order graph with these counter- 
terms has only integrable soft-gluon singularities since the singularities of the counterterms 
match the singularities of the cut order af +1 graph when either of the gluons become soft. 
Furthermore, the cancelling terms have the same final states: we are cancelling graphs with 
four final state partons against counterterms with the same four final state partons. 

B. Virtual soft gluons 

Subtracting the virtual terms, X[soft, virtual] requires a bit more analysis. We want to 
subtract these terms from the cut next-to-leading order graphs with three partons in the 
final state - a quark, a gluon, and an antiquark. Such a cut graph has a virtual loop 
either to the right of the final state cut or to the left. We want the contributions from 
— X[soft, virtual] to cancel the soft gluon singularities from these virtual loops. However, the 
terms in X[soft, virtual], as written, do not have the right structure to effect the cancellation. 
Thus we need to rearrange the terms. 

In the {i, j} term in X[so ft, virtual] we have the function 

-F l3 (i) = -c l3 (LMjk. (42) 

In the {j, i} term we choose to redefine the integration variable so that I — > —I. Then we 
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have a contribution with 



We use the identity 
where 



-F tJ (-t) = -C tJ ^MM . (43) 



-F^l) - Fy(-Z) = Fj<{l) + F«(l), (44) 



Qi ■ Qj - (i-Qi)(l-Qj) 



(Qj - Qi) - l + ie \l + qi-l l-qj-l, 



2Cij 



(45) 



and 



Qi ■ Qj - (l ■ Qi)(l • Qj) 



(Qi-Qj)-l-ic \l-qi-l 1 + qj-l, 



2C{j 



(Qi - Qj) - l-ie 



(46) 



The function gives the singular factor that arises from the exchange of a soft virtual 
gluon in a loop graph to the left of the final state cut. Similarly, gives the soft gluon 
singularity for a loop graph to the right of the final state cut. The essential part of a 
demonstration of this results from considering the integral 

j 2 .n3/g iCyu^D^i) 

1 1 J 2ir (u i -l + ie)(-u j -l + ie)(l 2 + ie)' 1 ' 

Here we have used the notation of Eq. (5) for the emission and absorption of a real gluon 
and have applied the same approximations. Upon performing the integration, we find 

I = Ft 3 il)- (48) 

In F?j(l), the two parts of the first term represent respectively the process in which parton 
j emits a (transversely polarized) soft gluon that is absorbed by parton i and the process in 
which parton % emits a soft gluon that is absorbed by parton j. The last term is the phase 
produced by the Coulomb force between the outgoing partons. 

With these functions at hand, we consider a cut next-to-leading order graph with three 
parton lines cut. Then there is a virtual loop. In the virtual loop there may be one or two 
gluons that connect two lines that go to the final state. Designate one of these gluons as 
potentially soft. This graph then corresponds to one of the contributions to X[soft, virtual]. 
We subtract a term of the form 



X-j [soft, virtual] = J dq\ J dq 2 Jdq 3 S(^2q^ 



^sof t d\i\ f d 2 i 

An 



X ^ F0; q u q 2 , q 3 ) R (q ± , q 2 , q 3 ), (49) 

7T J 
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with J = L or J = R depending on whether the virtual loop was to the left of the cut or to 
the right. We do this for each potentially soft gluon in the virtual loop. 

This procedure takes care of the terms in X[soft, virtual] that we needed to subtract from 
X. Additionally, the singularity of the counterterms matches the singularity of the graph 
when either of the gluons become soft (see [11]) so that the graph with these counterterms 
has only integrable soft-gluon singularities. 

There are some technical issues to attend to. First, one should match the integration 
variables in the loop graph to the integration variables in the counterterms, so that {q±, g*2, 93} 

— * 

are the momenta of the final state partons and I is the momentum of the potentially soft 
gluon in each case. Second, one needs to deform the integration contour for the counterterms 
according to the ie prescriptions indicated. This deformation should closely match that of 
the next-to-leading order graph in the limit of small /, as discussed in the Appendix. 

VI. SECONDARY SHOWERING FROM NLO GRAPHS 

We have started with the order af +1 graphs simply as specified by the Feynman rules 
and then, because there are order af +1 effects incorporated in the splittings of the partons 
emerging from the af graphs and in the soft gluon radiation from the af graphs, we have 
subtracted certain quantities from the af +1 graphs. The result is that the subtracted af +1 
graphs with and without virtual loops are separately free of infrared divergences. 

One can attach a shower to each final state parton emerging from a subtracted order 
af +l graph. Each shower affects the expectation value of the observable being calculated. 
However, according to the argument in Ref. [1], as long as the observable is infrared safe, the 
effect is to add terms proportional to af +2 and higher powers of a s . This holds independently 
of the exact functions used to generate the showers, as long as these functions have the 
required generic properties. For our present purposes, we can use the simple algorithm 
described in Sec. VIII of Ref. [1]. 

The only thing that we need to specify is the choice for width parameter k for the first 
splitting in the shower. We will take 



where q is the momentum of the parton that is to split and t is the thrust of the partonic 
final state in the order af +1 graph. The factor (1 — t) is important. In the case that we 
generate an event with thrust very near to 1, we do not want the showering to turn this 
two-jet event into a three-jet event, an event with thrust substantially less than 1. The 
factor (1 — t) prevents this. [For thrust greater than a value t max close to 1, the code [14] 
provides a stronger cutoff: k 2 = q 2 (1 — t) 2 /(l — t max ) 2 . The default value of t max is 0.95.] 

VII. EXPONENTIATION 

We now return to the terms that we have added. Consider the sum 



n 2 = q 2 (1-t), 



(50) 



J[Born] + Jfsoft, real] + J[soft, virtual] 
from Eqs. (22), (23), and (24). Compare this to the quantity 



(51) 
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X J2~ F ij( 1 ' 9l> 92, 5s) 9l, 52, 9s) 

'in i/ 7 



xexp Hn 77 ./ 77 ^ T ' 7; ' 72 ' ' < 52 > 



Here the three outgoing partons radiate a soft gluon according to the radiation pattern 

— * 

specified by the Fy. The gluon radiation at small I is suppressed by the exponential function, 
which gives the probability that the gluon was not already radiated with a higher energy than 
\l |. As we found in the companion paper [1] on collinear emissions, there is a subtraction to 
multiplication theorem that says that 

J[Born] + J[soft, real] + J[soft, virtual] = J[radiate] x (l + 0(a 2 )) . (53) 

The proof of this is simple. We have 

Xfradiate] = X[radiate, 1] + Xfradiate, 2], (54) 

where 

rMson d\l\ r d 2 l 



J[radiate, 1] = / dq\ J dq 2 J dq 3 5 (J2 Qi 



o \l\ J 4tt 



x^T F a(h 5i, 52, 9s) 5i, 52, qs) ~ Ro(qi, 52, &)} 

^soft d\i'\ r d 2 i 



x 



ex P ( - ^ / — £ - F,,,(/ ; 52, 53) I (55) 



and 



Jfradiate, 2] = Jdq\ J dq 2 J ' dq 3 S(j2vi) #0(51,52,53) 

x io "]ff 7 47 £ V y( ; 9i ' fe gs) (56) 

X exp \- ^ jjj—j:- F.Al ; q u q 2 , q 3 ) 
Now X[radiate, 1] can be expanded in powers of a s to give 

Jfradiate, 1] = J dq x J dq 2 J dq 3 5 (£ Qi) J q "jjy J 

x £ — 5i, 52, 53) [Rijik 5i, 52, 5s) - ^o(5t, 52, 53)] 

+£>(a 2 x (57) 

This is just X[soft, real] + X[soft, virtual]. The integral in Xfradiate, 2] looks complicated, but 
it is simple because the \l | integral is the integral of a derivative. Thus X[radiate, 2] is the 
difference of the integrand between the integration end points: 

X[radiate, 2] = f dq x J dq 2 J dq 3 S(j2 ft) R o(Qu 52, q 3 )[F(M soh ) - F(0)], (58) 
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where 



= exp (-(a a /7r) CTog(M soft /|f|)) . (59) 
Here C is the angular integral in the exponent, given in Eq. (25). Now F(M soft ) = exp(— 0) = 

— * 

I. On the other hand, the exponent diverges as |Z | — > 0. Furthermore, C is positive. Thus 
F(\l |) — ■> exp(— oo) = as \l \ — > 0. The result for X[radiate, 2] is then 

J[radiate, 2] = J dq-y J dq 2 J dq 3 S(j^ ft) #o(<fi, 92, %)■ (60) 
This completes the proof of the theorem. 



VIII. THE RESULT 

We are now able to generate a useful expression for the observable X. We begin with 
a cut graph that contributes to X[Born], Eq. (12). For each of the final state partons, we 
introduce an additional integration, 

representing the splitting of that parton, as explained in Ref. [1]. We also introduce an 
integration, 

f M so(t d\l\ f d 2 l a . A 

JO |/| ■/ 47T ^ 7T 

xexp (-/n IFf/tegT^*'*'*']' (62) 

representing the radiation of a soft gluon from the final state partons. The function repre- 
senting the Born matrix element then depends on all of the integration variables. This gives 
a probability to create a total of seven partons. Each of these seven partons then serves as 
the source of a secondary shower. 5 This is illustrated in Fig. 2. 

In addition, there are contributions from order af +1 cut Feynman graphs plus their soft 
gluon counterterms. Included here are cut self-energy graphs, but only for virtuality above 
the cutoff Aygf , so that there are no collinear divergences. Each of the three of four partons 
emerging from such a graph then serves as the source of a secondary shower. These order 
af +1 remnant terms are complicated, but they have only integrable singularities and they 
are just what is needed to leave the calculation correct to next-to-leading order. 



5 For the soft gluon, we choose the starting scale for the secondary shower to be k 2 — c KlSO ft / 2 where I is 
the momentum of the soft gluon and c K SO f t is a constant taken to be 1/10 in the numerical example in 
the following section. 
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IX. A NUMERICAL TEST 



In this section I present some numerical results based on the algorithm described above. 
I have computed the thrust distribution da/dt for thrust t = 0.86, in the middle of the 
three-jet region. We will compare this distribution calculated with parton showers to the 
same distribution calculated with a straightforward NLO computation with no showers. The 
difference should be of order af +2 - We divide the difference by the NLO result, forming 

(NLQ-shower) - NLO 
R = NLO • (63) 

The ratio R should have a perturbative expansion that begins at order a 2 . We can test this 
by plotting the ratio against a 2 . We expect to see a curve that approximates a straight line 
through zero for small a 2 . For comparison, I exhibit also the ratio with the order graphs 
with showers but with the af +1 corrections omitted, 

(LO-shower) - NLO 

= nl5 ' (64) 

The ratio .Rlo should have a perturbative expansion that begins at order a]. Thus we expect 
to see a curve proportional to the square root function [a 2 ] 1 / 2 for small a 2 . The size of the 
coefficient of [a^] 1 ^ in R LO does not have any great significance, since it is quite sensitive 
to the choice of renormalization scale. 

The results of this test are shown in Fig. 3 for = Mz with a renormalization scale 
fi = y/S/6. The values of a s {M z ) range from (1/2) x 0.118 to 2 x 0.118, where 0.118 
represents something close to the physical value for a s (Mz)- One should note that 2 x 0.118 
amounts to quite a large a s in this calculation since a s (M z ) = 0.236 gives a s (fj) = 0.586. 

We see the expected shape of the -Rlo curve. We then examine whether the R curve 
approaches a straight line through the origin as a 2 — > 0. Within the errors, it does. However, 
the slope of the straight line is quite small. Presumably with other choices of parameters in 
the program, the absolute value of the slope would be bigger. 

In Fig. 4, I show the same comparison, but this time for t = 0.71. This is near the 
value t = 2/3 that marks the far end of the three-jet region, with the three partons in what 
is sometimes called the Mercedes configuration. The results are similar to the results for 
t = 0.86. 

In Fig. 5, I show the same comparison, but this time for t = 0.95. This is near the two-jet 
limit at t — 1. For a s = 0.118, one would normally not use a calculation that did not include 
a summation of logs of 1 — t for t this close to 1, since log(0.05) 2 ~ 9. Thus I would not 
recommend using the code discussed in this paper for a comparison to data this near to the 
two-jet limit. Nevertheless, we can still test for the absence of an a I term in R. Looking at 
the graph, we see that, within the errors, there is no evidence for a nonzero a] term in R. 
We expect an a 2 term, but it appears that the coefficient of a 2 is quite small. On the other 
hand, it appears that some of the yet higher order terms are quite substantial. 



X. AREAS WHERE MORE WORK IS NEEDED 

The present work suffers from a number of deficiencies. First, the code is designed to 
calculate three-jet cross sections in electron-positron annihilation correctly to NLO, but it 
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FIG. 3: Comparison of the NLO calculation with showers added as described in this paper and 
[1] to a pure NLO calculation using [14]. I plot the ratio R defined in Eq. (63) for the thrust 
distribution at thrust equal 0.86. Also shown is the ratio -Rlo ; defined in Eq. (64), in which the 
order af +1 correction terms are omitted from the calculation. The cm. energy is = Mz and 
the renormalization scale is chosen to be /j = y/S/G. These ratios are calculated for a s (Mz) 2 = 
{0.25,1,2,3,4} x (0.118) 2 and plotted versus a s (M z ) 2 /(0.118) 2 . 
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FIG. 4: Comparison of the NLO calculation with showers added to a pure NLO calculation for 
t = 0.71. The notation is as in Fig. 3. 
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FIG. 5: Comparison of the NLO calculation with showers added to a pure NLO calculation for 
t = 0.95. The notation is as in Fig. 3. 

ignores the two-jet cross section. It should be rather simple to calculate the two-jet cross 
section at next-to-leading order with showers. Then one should merge the two calculations 
so that three-jet observables are calculated correctly at order a 2 while the two-jet cross 
section is calculated correctly at order a\. 

Such a merged calculation would have benefits even for the calculation of three-jet ob- 
servables. Specifically, there would be an advantage when calculating three-jet observables 
close to the two-jet region. For example, consider the thrust distribution near to t — 1. 
The two-jet calculation with showers, before any merging, would approximately sum terms 
proportional to log J (l — t), generating these terms by showering from the two-jet hard 
process. For da/dt, the merged calculation should contain the a] and a 2 contributions cor- 
rectly, including of course the term proportional to a 2 log(l — t). It should also include the 
terms a™ x logs for n > 2 generated by showering from a two-jet configuration. The three-jet 
algorithm presented in this paper does not include merging with a two-jet calculation and 
thus does not generate these logarithms. Thus the present program, when applied to the 
thrust distribution, should not be used for small values of 1 — t. 

Second, it is desirable to add not only parton showers but hadronization to NLO per- 
turbative calculations. For hadronization, one can use an existing model [5-7]. However, 
most hadronization models require color information on the final state partons. Since the 
current code does not assign a color structure to these partons, it will be necessary to add 
this capability. 

Third, in principle, hadronization models should be infrared safe in the sense that final 
partonic states that differ by a parton splitting into two almost collinear partons (or one 
parton and another with almost zero momentum) will produce almost the same hadronic 
final state. Since hadronization models currently in use do not necessarily have this property, 
the NLO calculation should provide it by combining parton pairs with virtuality less than 
some cutoff parameter matched to the hadronization model to be used. 

Finally, the current paper includes hardly any numerical investigations to check how well 





$ \ 






i • 




■ 

f B R LO 




[£=0.95 i 


1 



20 



the program works for various observables and how the results depend on parameters such 
as the renormalization scale /i and the soft gluon limiting energy M so f t . Addressing these 
needs remains a topic for future research. 

XI. CONCLUSIONS 

In this paper and [1], we have seen how to add parton showers to a next-to- leading order 
calculation in QCD in such a way that the result obtained for an infrared safe observable 
remains correct to next-to-leading order. The main feature of this procedure is to turn the 
collinear singularities of the order af +1 graphs into the primary parton splittings, that is the 
first splittings of the partons that emerge from an af graph. This was accomplished in [1] in a 
particularly simple fashion by virtue of the fact that Ref. [1] is based on a NLO calculation in 
the Coulomb gauge. In this gauge the collinear and collinear x soft singularities are entirely 
contained in the cut self-energy graphs. Thus we simply eliminate the small-virtuality part 
of cut self-energy subgraphs in af +1 graphs and instead incorporate these contributions 
into the primary parton splittings from the corresponding af graph. One could subtract a 
different quantity from the af +1 graphs as long as it has the same collinear and collinear x 
soft singularities, 6 as in the work of Refs. [3, 4]. 

The next step, and the main subject of this paper, is to treat the soft, wide angle sin- 
gularities of the theory. Again, we subtract terms with these singularities from the af +1 
graphs and then incorporate them as radiation from the af graphs. Here, the hard partons 
emerging from an af scattering form an antenna that radiates the soft gluon. The gluon is 
radiated coherently from the antenna and is not part of any jet. 

This gives an algorithm that can be viewed as a simulation of the production of three 
hard partons (for the three-jet process considered in this paper) and one soft gluon, with 
each of the three hard partons splitting into two daughter partons. The six daughter partons 
and one soft gluon then serve as the progenitors of parton showers according to the simple 
algorithm used in Ref. [1] or any standard algorithm. There remains the remnants of the 
order af +1 graphs, with three and four parton final states. Each of these final state partons 
serves as the progenitor of whatever sort of parton shower is desired. 

The results of this paper suggest a strategy for building future NLO calculations matched 
to parton showers. The matching concerns the primary parton splittings, that is the first 
splittings of the partons emerging from a Born graph. The matching also needs the radiation 
of a soft gluon from the antenna formed by partons emerging from the Born graph. This 
suggests that the author of an NLO calculation might build in the primary parton splittings 
and the soft gluon radiation. There are also secondary splittings. These are the splittings 
from the soft gluon produced from a Born graph, the splittings of the daughters produced by 
the primary splittings, the splittings of the partons from an order af +1 graph, and, finally, 
the further splittings of all of these partons as a complete parton shower is formed. The 
secondary splittings do not need to be matched to the NLO calculation. That is, the soft 
gluon and the partons after the primary splittings and the partons produced by the af +1 
graphs can be fed to any shower and hadronization algorithm as long as it is infrared safe. 
The result is that for a given process such as e + + e~ — > 3 jets we can have ./Vnlo next-to- 



This might, however, require adjustment of the treatment of the soft gluon singularities presented in this 
paper. 
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leading order calculations and Nmc shower and hadronization calculations and we do not 
need Nnlo x ^mc matching schemes as long as each NLO calculation includes its primary 
parton splittings and treatment of soft gluon radiation. 

This does not preclude debate over how best to do the primary parton splitting and soft 
gluon radiation. This paper and Ref. [1] have, of necessity, chosen one scheme, but that 
scheme may well be less than optimal. 

The code described in this paper is available at [14]. 
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APPENDIX: THE INTEGRATION FOR VIRTUAL SOFT GLUONS 

We have seen that in order to calculate the contribution to X from an order a 2 s graph 
with a virtual loop, we should calculate a difference integral of the form 

J dqi J dq 2 J dq 3 {J d 3 l f(l, q u q 2 , <f 3 ) 

-E / ^Lw* < Ml it )^F0' ]qi ,q 2 ,q 3 )R o (q u q\q 3 )\, (A.l) 

where J is L or R and F£j and F$ are given in Eqs. (45) and (46). In this appendix, we 
consider J = L, that is a virtual loop to the left of the final state cut. The integrand for the 
original graph is represented here by the function /. The momenta are the momenta of 
the final state partons, while / is the momentum in the virtual loop. There are either one 
or two terms in the sum over the indices {i,j}, depending on the graph. The term {i,j} 
is included if there is a potentially soft virtual gluon connecting final state lines i and j in 

— * — * 

/. Let lij be that combination of I and the q^ that is carried by this potentially soft virtual 
line in /. For example if the final state propagators 1 and 2 are connected by a gluon line 
that carries momentum I in the virtual loop and propagators 2 and 3 are also connected 
by a virtual gluon line that carries momentum I — q 2 , then l 12 = I and l 23 = I — q 2 . In the 

— * — * — * 

virtual subtraction term, we call the integration variable I'. If we identify I' with 1^, the F^ 

— * 

subtraction term matches / when 1^ — > 0, so that the leading singularity is cancelled. 

The method used in this paper is to perform the virtual loop integrations numerically, 
as in [2, 10-12] and [1]. Thus we must make sure that this cancellation works numerically. 
First, we define the dependence of / and R on qi,q 2 ,q 3 to be the dependence computed 
from the relevant Feynman graphs with the final state momenta qi,q 2 ,q 3 on-shell: q\ = 0. 
Of course, in our present discussion, the qu are the momenta of on-shell massless partons, 
so evidently q\ — 0. However, at a later stage in the calculation each of these partons will 
generate a parton shower with the same q. Even with a shower, we take the functions / and 
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R to remain unchanged. This is in contrast to our treatment of showering from the Born 
graphs and is required for the NLO graphs in order to preserve the soft gluon cancellations. 

Next, in the integral of the full graph /, the integration is deformed 7 into complex loop- 
momentum space, 

/ d 3 lj(l) f(l + m(f), gi, q 2 , q 3 ) . (A.2) 

The imaginary part of the deformed loop momentum is a definite function k(1) of the real 
part. The Jacobian J is the determinant of the matrix (5u + idn 1 /dl J ). By deforming the 
contour, we avoid the singularity of the form 

V => , (A.3) 

\Qi\ + \Qj\ ~ W + hj + i*\ - \Qj - kj - in\ + it 

in /, the {i,j} scattering singularity. It is, however, not allowed to deform the contour away 
from the soft singularity at kj = 0. For that reason k(1 ) is proportional to /|- as ff- — > 0. 

Hence we come very near to the scattering singularity when is small. 
What about the subtraction term? Here there is a singularity of the form 

^ . (A.4) 

V ■ [Qj ~ Qi) + ^ 

We may call this the approximate scattering singularity. We note that 

Iftl + \Qj\ ~ Wi + l v\ ~ \<U ~ l v\ ~ Uj ■ (Qj - Qi) (A. 5) 

— # — # 

for 1^ — > 0, so the singularity (A.4) is indeed an approximation to Eq. (A.3) if we identify /' 
with lij. As in the / integral, we need to deform the contour to avoid the singularity (A.4). 
As in the / integral, the deformation will have to vanish as we approach the soft singularity 
at I' — 0. In order for the cancellation between the full graph and the soft gluon subtraction 
to work, we need to be careful about the matching of the singularities in the two integrands. 

We are now in a position to state the matching problem more precisely. We define complex 
vectors s^-, one for each potentially soft gluon line {i,j} in Eq. (A.l). Each s^- is a function 
of the corresponding 1^ that is asymptotically equal to 1^ when 1^ is small. Then we write 
the integral in Eq. (A.l) in the form 

Jdqi Jalq 2 Jdq 3 5(j2q^ J d 3 l ^J(l) f(l + i«, cfi, q 2 , <f 3 ) 

" E f|4 < <ft) ^ l r fMr ft, ft, &) & &)}■ (A.6) 

The Jacobian J s is the determinant of the matrix dsfj/dl^. Our goal will be to specify the 
functions so that we avoid the approximate scattering singularity in the second term of 
Eq. (A.6) and so that the 1^ — > singularities in the sum of the two terms cancel sufficiently 
well that the integral is convergent. 

— * 

We will approach this problem in stages. First, we concentrate on the 1^ — > region. 



7 Some relevant theorems concerning contour deformations in more than one dimension are given in [11]. 
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We immediately recognize a small problem. The soft singularity (with k — 0) lies on an 
ellipse, while the approximate soft singularity lies on a plane. To make them match at the 
required accuracy, we define 

% = kj + (1 - Z~J/M s 2 oft ) {£(4) + iCikj) } (to be improved), (A.7) 

where £ and £ are certain functions of kj to be defined below and M s 2 oft is the maximum 
value of / 2 in the soft gluon subtraction, Eq. (14). 

The idea is to have the soft subtraction cancel the leading singular behavior of / point 

— * 

by point in I. To do this, we need to match the singularity 

(A.8) 



Do(hj + |ft| + - |ft + kj + i%\ - Wj - kj + i£ 

in / with the singularity 



D a (sij) Sij-iqj-qi) 



' (A.9) 



in Fij. 



There are three parts in the relation (A.7) between Sij and kj- There is a real displacement 
£, a contour deformation ( in the imaginary direction, and a factor (1 — /j 2 /M 2 oft ). The factor 

(1 — / 2 /M 2 oft ) serves to set the deformation to zero at the edge of the integration region and 
needs no further discussion. We choose the real displacement to be 

lift) = oLriLrirf* ,i {t 1 " ('« • ft) 2 ] ft " t 1 " & ' ^) 2 ] ft} • ( A - 10 ) 

z lftl |yj| L 1 % ' yjj 

This choice is motivated by carrying out the expansion (A. 4) to one more order. Then one 
finds that 

(kj + i(kj)) ■ (Qj ~ ft) ~ Wi\ + Wj\ - | ft + 4'l - Ift - kj\ + °^j)- ( A - n ) 
The imaginary displacement is 

<f(f) = Cc f 2 min((l + Z • (1 - Z • qj)) (q, - q t ). (A.12) 

Here is a parameter with dimensions of inverse momentum. This function matches the 
function k(1 ) in the limit of small I as long as one chooses to be the same parameter as 
is used in k(1). s Thus 

DMj) = D (kj + »«) + 0(1^). ( A - 13 ) 
There are other singularities that have to match. First, both / and the soft subtraction 
terms have l/l^l 3 soft singularities that evidently match since % = kj to lowest order. 
Second, both terms have collinear singularities along the lines l+Z^-^j = and 1 — kj-qj = 0. 
These are pinch singularities: it is not allowed to deform the contours to avoid them. For 



Specifically, — 2a/[(l +7)(|%| + \qj\ — \q% + <fj|)], where a and 7 are the parameters used to define the 
deformation in [11]. 
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this reason ((I) vanishes when I is in the same direction as — <fj or qj. Furthermore, the part 

— # — * 

of £ that is not proportional to qi vanishes when / is in the same direction as — qi and the 
part of £ that is not proportional to g} vanishes when I is in the same direction as cfj. 

We are now ready for the second stage of our construction. In the first stage, we con- 

— * 

centrated on the ly — > region. The function s^- given in Eq. (A. 7) works fine for small 

— * — * 

lij. However, our integration region ends at = M so f t , so we supplied a factor that turns 
the contour deformation off at the boundary of the integration region. 9 Unfortunately, the 

— * — * 

singularity at kj ■ (qj — qi) =0 intersects the surface \kj\ = M so f t . Thus we should retain 
some contour deformation even at the boundary \kj\ = M so f t of the integration region. This 
is simple to do. 

Let us consider the problem anew, this time ignoring the special problems at \l | — > that 

— * 

were addressed in the analysis above. Let n = (cjj — qi)/\qj — q%\ and write kj in spherical 
polar coordinates (I, 9, <fi) based on n as the z axis. Then we have an integral 

/ l 2 dl dx d<j), (A. 14) 

JO J-l J-TT 

where x = cos 9, that is 

x = kj ■ n/\hj\. (A.15) 

The integrand has a singularity of the form l/(x + ie). All that we need to do is deform 
the contour for the x-integral away from the singularity. There are collinear singularities at 
lij ■ (jj = 1 and kj ■ cji = —1, that is at <fi = 0, n and x = x , where 



x = n-q j = y/(l-q i -q j )/2. (A.16) 
The contour cannot be moved at the collinear singularities. Thus we can take 

x — > z = x + iy, (A. 17) 

where 

V = (?/M s 2 oft ) (x 2 - x 2 ) 2 9(x 2 < xl). (A.18) 

The factor /?/M s 2 oft is included here simply to turn this deformation off for small Ifp since 
we already have a method for dealing with the small Ifj region. In a coordinate frame 

— * 

independent notation, this is 1^ — > % with 

% = kj + v(hj) (to be improved), (A. 19) 

where 

ff=iy \lij\n+ I J |_^ 2 - 1 J (hj ~ (kj • n) nj . (A.20) 



This is the analogue of keeping the contour fixed at the end point of the integration in a one dimensional 
contour integration. 
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Now we have two methods for moving the integration contour, each of which works in its 
region and turns off in the region of the other method. All that we have to do is to add the 
two deformations, 

% = kj + (i - 5/M s 2 oft ) {ify + <(/;.)} + rffy. (A.21) 

This is the final prescription, with £ given in Eq. (A. 10), ( given in Eq. (A. 12), and fj given 
in Eq. (A.20). 
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